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We present an accurate method of treating the one-neutron removal reaction at interme- 
diate incident energies induced by both nuclear and Coulomb interactions. In the method, the 
nuclear and Coulomb breakup processes are consistently treated by the method of continuum 
discretized coupled channels without making the adiabatic approximation to the Coulomb 
. interaction, so that the removal cross section calculated never diverges. This method is 

' applied to recently measured one-neutron removal cross section for ^^Ne+^'^C scattering at 

^ I 230 MeV/nucleon and ^^Ne-|-^°*Pb scattering at 234 MeV/nucleon. The spectroscopic factor 

and the asymptotic normalization coefficient of the last neutron in "^^Ne are evaluated. 

§1. Introduction 

. Unstable nuclei have exotic properties such as the halo structur^^SJ and the 

change of magicity for nuclei in the region called "Island of inversion" S"® These 
novel quantum properties have inspired extensive experimental and theoretical stud- 
ies. An important experimental tool of exploring such exotic properties is the nucleon 
\ removal reaction; see for example Ref. The theoretical tool of analyzing such an 

^ ■ inclusive reaction is the Glauber modelP The theoretical foundation of the model is 

investigated in Ref. Ej) . The model is based on the eikonal and the adiabatic approx- 
imation. It is well known that the adiabatic approximation makes the removal cross 
section diverge when the Coulomb interaction is included. The Glauber model has 
\(Q . thus been applied mainly for lighter targets in which the Coulomb interaction is neg- 

\ ligible; see for example Refs. [7|) .[T0 |) -[15 ]) . Very recently, the Glauber model with the 

' Coulomb correction was proposed.'^''^ In the Coulomb-corrected eikonal model 

(CCE),II^'II3 the El (dipole) component of the eikonal Coulomb phase is replaced 
(<~^ \ by that based on the first-order perturbation. For the elastic-breakup scattering of 

^^Be from a 208p|j target at 69 MeV/nucleon, it was shown in Ref. [I?]) that CCE 
well simulates results of the dynamical eikonal approximation (DEA).'^ 

Meanwhile, the method of continuum discretized coupled channels (CDCC^(^''23 
^ ■ is an accurate method of treating exclusive reactions such as the elastic scattering and 

, the elastic breakup reaction in which the target is not excited whereas the projectile 

is broken up into its fragments. The theoretical foundation of CDCC was shown in 
Refs. [2T ]) -l23 ]l . Actually, CDCC has succeeded in reproducing data on the scattering 
of not only stable nuclei but also unstable nuclei; see for example Refs. [M|) -[30 |) and 
references therein. The dynamical eikonal approximation (DEA)'^ is also an accu- 
rate method of treating exclusive reactions at intermediate and high incident energies 
where the eikonal approximation is reliable. The nucleon removal reaction is com- 
posed of the exclusive elastic-breakup component and the inclusive nucleon-stripping 
component. CDCC and DEA can evaluate the elastic-breakup cross section, but not 
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the stripping cross section. 

The experimental exploration of halo nuclei is moving from lighter nuclei such 
as He and C isotopes to relatively heavier nuclei such as Ne isotopes. Very recently, 
a halo structure of ^^Ne has been reported by the experiment on the one-neutron 
removal reaction a-n at 230 MeV/nucleon not only for a ^^C target but also for 
a ^°^Pb target.!^ This is the heaviest halo nucleus in the present stage confirmed 
experimentally and also reside within the region of "Island of inversion" . The deter- 
mination of the spin-parity of ^^Ne in its ground state is essential to understand the 
nature of "Island of inversion" . The one-neutron removal reaction has been analyzed 
with the Glauber modelpSt for a ^'^'^Pb target, the elastic breakup component due 
to one-step El transition is added. 

In this Letter, we present an accurate method of treating the one-neutron re- 
moval reaction at intermediate incident energies such as 200 MeV/nucleon induced 
by both nuclear and Coulomb interactions. In the method, the nuclear and Coulomb 
breakup processes are accurately treated by CDCC without making the adiabatic ap- 
proximation to the latter, so that the removal cross section calculated never diverges 
and reliable even in the presence of the Coulomb interaction. Thus, this method is an 
essential extension of the Glauber theory and CDCC. This method is applied to the 
one- neutron removal reaction of the ^^Ne-I-^'^C scattering at 230 MeV/nucleon and 
the ^^Ne-|-^^^Pb scattering at 234 MeV/nucleon. The results of the present method 
are compared with those of CCE. The spectroscopic factor S and the asymptotic 
normalization coefficient Canc'^'^'* of the last neutron in '^^Ne are evaluated with the 
new theory. 



§2. Formulation 



We consider the one-neutron halo nucleus as a projectile (P) and take a single- 
particle model for the nucleus; namely, the projectile consists of a core nucleus (c) 
and a neutron (n). The scattering of P on a target (T) is then described by the 
three-body (c-|-n-|-T) Schrodinger equation 



— V| + h + U{r, 



E 



ip- = 



il) 



with the interaction 



C/=[7f)(r„) + [7f)(r,) + [/f)(., 



where fi is the reduced mass between P and T. The three-dimensional vector R = 
(6, Z) stands for the coordinate between P and T, whereas (x=n or c) is the 
coordinate between x and A and r means the coordinate between c and n. The 
operator h = Tr + y{r) is the projectile Hamiltonian composed of the kinetic-energy 

operator and the interaction V{r). The potential is the nuclear part of the 

(c) 

optical potential between x and T, whereas Uc is the Coulomb interaction between 
c and T. 

Now, we consider scattering at intermediate incident energies, say 200 MeV/nucleon, 
for which the eikonal approximation is quite good. The three-body wave function 1^ 
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is first assumed to be 



= 0^{R,r), 



(3) 



where the operator O is defined by 



O 



JK-Z 



(4) 



with the wave-number operator K = y^2ij,{E — h) /h and the velocity operator v = 
hK / of the relative motion between P and T. When Q is inserted into ([T]), we have 
three terms 0[V^^], [V^jO] • and [V|jO]V'- The first term is much smaller 

than the others, since tp is slowly varying with R compared with O. Neglecting the 
first term leads to 



aZ 



(5) 



Regarding Z as "time" virtually and solving ^ iteratively, we obtain the formal 
solution 



■0 = exp 



iV / dZ'O^UO 



(6) 



where V is the "time" ordering operator. Taking Z to oo in ([6]), we finally get the 
S'-matrix operator 



S = exp 



iV / dZO^UO 



(7) 



In the Glauber model, the adiabatic approximation is made as the secondary ap- 
proximation. In the approximation, h is replaced by the ground-state energy eo, and 
hence O^UO and V in 1^ are reduced to U/{hvo) and 1, respectively, where vq is 
the velocity of P in the ground state relative to T. This is nothing but the S'-matrix 
in the Glauber model. 

At intermediate incident energies, the adiabatic approximation is good for the 

but not for the long-range Coulomb 



short-range nuclear interactions, u2^^ and Uc^^ 



(c) 

interaction U} ■ This can be understood from the matrix element 



dz{ip^\6^u6\ip^) 



JiKo-K)Ru foo 



hvo 



dZ{ipf.\U\ipo) 



(8) 



between the ground state tpo of P with the intrinsic energy eo and the continuum 
state (pf^ of P with the intrinsic momentum and energy, hk and e{k). In ([8|), HKq 
(hK) is the momentum of P in the ground (continuum) state relative to T, and 
Rjj is the range of the interaction considered. For the ^^Ne breakup reaction at 
200 MeV/nucleon, the typical excitation energy e{k) ~ 1 MeV. For the Coulomb 
interaction, because of its long-range property, {Kq — K)Ru is large even if Kq — K 
is small. This means that we can not set Kq — K lo zero, that is, the adiabatic 
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approximation does not work. Actually, the elastic-breakup cross section diverges in 
the adiabatic limit of K = Kq because of the slow decrease of the matrix element 
([S]) in b^^^^ For the nuclear interactions, meanwhile, Ru ~ 11 fm and hence 
{Kq — K)Rij ~ 0.074. Thus, the adiabatic approximation is good for the nuclear 
interactions. This indicates that U^'' is commutable with O with high accuracy. 
Therefore, we can take the replacement 

otf/f)d^c/f)/(/.^o). (9) 

The replacement ^ is accurate as shown later by numerical calculations. Using this 
replacement, one can get the important result 



S — SiiSc 



with 



5n = exp 



hvQ 



X) 

Sc = exp -iV I dZdt(c/jN) + uP)d . 

J —oo 

Equation ([TO]) can be derived as follow. The replacement ([9]) leads ([5]) to 



.dtp 



with Uc 



uP + uP. 



Defining as ip 
Q = exp 



with 



hvo 



and inserting ip 



into (|T3|) leads to 



(10) 

(11) 
(12) 

(13) 

(14) 
(15) 



where use has been made of the fact that Q, i.e. uP , is commutable with O with 
high accuracy. The formal solution of p5p is 5c, whereas is obtained from Q by 
taking Z — )• oo. Noting that ^ = Q^, one can reach (fTUj) . 

Thus, S can be separated into the neutron part and the core part ^c. We 
can not calculate directly with (jl2p . because it includes the operators O and V. 
However, 5c is the approximate solution of the Schrodinger equation 



'2/x 



l^n = 0, 



(16) 



when the eikonal approximation is made. We can calculate the matrix elements 
((/?o|5c|9?o) and (99^ 1 5c I (/Jo) by solving ([16]) with eikonal-CDCC'^^J in which the eikonal 
approximation is made in the framework of CDCC. Non-eikonal corrections to 5c 
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can be easily made by solving ()16p with CDCC instead of eikonal-CDCC, although 
it is not necessary for the present intermediate scattering. As mentioned above, Sn 
is obtained by pT|) . 

The reaction theory constructed above is referred to as the eikonal reaction the- 
ory (ERT) in this Letter. We can test the accuracy of the eikonal approximation 
by comparing the S'-matrix elements of CDCC and eikonal-CDCC.'® We confirmed 
that the approximation is quite good for the present scattering system. This approx- 
imation is good even for lower incident energies such as 30 MeV/nucleon.® So all 
CDCC calculations in this Letter are done with eikonal-CDCC. 

We can derive several kinds of cross sections with the product form (jlOp . fol- 
lowing the formulation on the cross sections in the Glauber model.'^''^ The one- 
neutron removal cross section o"_„ is the sum of the total elastic breakup (diffractive) 
cross section cjbu and the total neutron-stripping cross section cjstr in which n is ab- 
sorbed by T: 

CT-n = CTbu + fstr (17) 

with 

0"str = J (fb{ipo\\Sc\'^{l - \Sn\'^)\<po) 

= [o-R - 0-bu] - [o-R(-n) - CJbu(-ra)], (18) 

where ctr and o"bu are the total reaction and elastic-breakup cross sections, respec- 
tively, defined by 

aK = j d^b[l - \{^o\ScSM\% (19) 
f^bu = / d^b[{^o\\ScSnf\^o) - \{^o\ScSM\% (20) 



and crR(— n) and cTbu(— "ra) correspond to the total reaction and elastic-breakup cross 
sections, respectively, in which ScSn is replaced by Sc- The last form of (fTSj) means 
that (T_„ can be obtained from cjr, cJbu) crR(— n) and crbu(— ''^) by solving ([1]) and 
(USD with CDCC. 



§3. Model setting 



We take the same three-body model as in Ref. [32]) except for minor difference, 
since the Glauber model calculation based on the three-body model yields a reason- 
able value of cr_n. The nuclear potential uji^^ (x=n or c) is calculated by the folding 
model in which the effective nucleon-nucleon (NN) interaction t is folded with the 
densities of particles T and x. As for t, we take the parameter set in Refs.[32]). l35p for 
E'nn = 240 MeV, where E'nn is the energy of NN scattering 0- The folding potential 
reproduces the reaction cross section of ^^C + ^^C scattering in the wide range of 

*' In Ref. l35|) . the efTective interaction is presented as the profile function _r(6NN)- It is a function 
of not the distance tnn between two nucleons but the impact parameter &nn- The interaction f(rNN) 
is so constructed that the profile function obtained from t(rNN) can agree with that of Ref. .35) . 
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incident energies!^ In Ref. [32]), the ^'^Ne density is calculated by the single-particle 
model in which the separation energies of last neutron and proton reproduce the ex- 
perimental values 3 MeV and 2.4 MeV, respectively; for the potential parameters, see 
the case of (rg, a) = (1.25, 0.75) in Table I of Ref. '32p, where tq and a are the radius 
and the diffuseness parameters of the Woods-Saxon form in units of fm. The target 
density is evaluated by the Hartree-Fock calculation with the effective NN interac- 
tion Gogny-DlS.'2S'E3 The ^'^Ne-n potential V is determined so that the neutron 
separation energy becomes 0.33 MeV; we take the case of (ro,a) = (1.25,0.75) in 
Table II of Ref. i32j) for lp3/2 and the case of (ro, a) = (1.25, 0.70) for 0/7/2. A large 
o"!^ indicates that ^^Ne is a one-neutron halo nucleus. This structure is realized 
with small angular momenta between c and n. However, the possibility of the 2sl/2 
orbit is small from the theoretical point of view, since the single-particle energy of 
2sl/2 is much higher than those of lp3/2 and 0/7/2. Hence, we consider the lp3/2 
and 0/7/2 orbits. 

§4. Results 

In ERT, the adiabatic approximation is assumed to be good for the nuclear po- 
tential Ui^\ This can be tested by setting U = ui^\rn) + uP (R) + uP{R) in 
the Schrodinger equation ([1]). In this setup, the projectile breakup is induced only 

by uP{rn), since the argument rc of uP and uP is replaced by R. Switching 
the adiabatic approximation on in the Schrodinger equation corresponds to the re- 
placement For the scattering of ^^Ne(lp3/2) from ^"^^Pb, the error due to the 
approximation is 0.2% for ctr, 1.9% for cjbui 4.1% for Ustr and 3.3% for (j_n- Thus, 
the error is small. 

Table |T] presents several kinds of cross sections and the spectroscopic factor S = 
^exp^^th^ for and ^o^Pb targets. As an important result, cS[lp3/2] =0.693 (0.682) 
for (208pb), whereas 5[0/7/2] =2.47 (5.65) for (208pb). Thus, S[lp3/2] little 
depends on the target and less than 1, but 5[0/7/2] does not satisfy these conditions. 
In Ref. I3ip . the Coulomb component of o"_n[2'^^Pb] for a ^osp]^ target is estimated 
to be 540 mb from the experimental values of cx-nP^^Pb] and 0"_n[i2C]. In ERT, the 
Coulomb component of (T_nP'^*Pb] agrees with o"bu[^'^^Pb] with good accuracy. The 
spectroscopic factor evaluated from the Coulomb component is S' = 5A0/a\^ = 0.675 
for the lp3/2 orbit and 7.36 for the 0/7/2 orbit. Thus, S' is consistent with S only 
for the lp3/2 orbit. Hence, we can infer that the major component of the ^^Neg.s, 
wave function is ^'^Ne(0"'") lp3/2 (5 ~ 0.69). We adopt in the following this 
configuration. 

In Ref. [32]) . S is calculated with the Glauber model; for a ^ospj-, target, the 
elastic breakup component due to one-step El transition is added. The resulting 
spectroscopic factor for the lp3/2 orbit is 5 = 0.822 for a target and 0.624 for a 
^''^Pb target The difference between the results of the present study and Ref. [32|) 
for mainly comes from the neglect of cxbu in the latter, and that for ^'^^Pb seems 
to stem from Coulomb higher-order contributions to cibu- Neglect of the Coulomb 
interaction in cistr may also be responsible for the difference. Thus, CCE is fairly 
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Table I. Several kinds of cross sections and the spectroscopic factors for Ne+^ C scattering at 
230 MeV/nucleon and ^^Ne+^°**Pb scattering at 234 MeV/nuclcon. The cross sections are 
presented in units of mb and the data are taken from Ref. I31|) . 
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90 


29 




244 


53 




CT-n 


114 


32 


79 


1044 


126 


712 


5 


0.693 


2.47 




0.682 


5.65 





good also for the present system. 

The potential V between c and n is not well known. Hence, S has a theoretical 
error coming from the potential ambiguity. The error is often estimated by changing 
each of rg and a by 30%. When the one-neutron separation energy of ^^Ne is 
0.33 MeV, S = 0.693 ± 0.133 ± 0.061 for a target and 0.682 ± 0.133 ± 0.062 for 
a 208p]2 target, where the second and third numbers following the mean value stand 
for the theoretical and experimental uncertainties, respectively. Thus, S includes 
a sizable theoretical error. This situation completely changes if we look at the 
asymptotic normalization coefficient (ANC). For r out of the range of V, ANC Canc 
is defined 

Iij{r) = CANchl^Hif^r) (21) 

with the radial part of the overlap function Iij{r) = {(j)c4'n\<P'p) , where (pc and (p^-^ are 

the intrinsic states of c and n, respectively, whereas 4>^ is the ground state of P. 
is a spherical Hankel function and n is the relative wave number between c and n in 
the ground state of P. For large r, Iij{r) is related to the normalized single-particle 
wave function ^pij{r) that is determined from V as 

Ii, (r) = VSipij (r) = ^/5CS^C^!+) (iKr) (22) 

with the sing le-particle ANC cjj^}.. Hence, we have 

Canc = V^cJ^;,. (23) 

When = 0.33 MeV, Canc = 0.320 ± 0.010 ± 0.028 fm-^/^ for a target and 
0.318 ± 0.008 ± 0.029 fm"^/^ fo^. 208^-^^ target. Thus, Canc has a much smaller 
theoretical error than S. This means that the one-nucleon removal reaction is quite 
peripheral. 

The experimental value of is not precisely determined: = 0.29ibl.64 MeV-l^^ 
It is thus better to see B^ dependence of Canc and S. When B^ = 0.1 MeV, 
Canc = 0.128 ±0.003 ±0.011 fm-^/^ and S = 0.530 ±0.084 ±0.047 for a target, 
and Canc = 0.105 ± 0.004 ± 0.010 fm-V2 and S = 0.358 ± 0.057 ± 0.033 for a ^ospb 
target. These values are plotted in Fig. [H Canc and S are sensitive to the value 



8 



Letters 



of Byi. We can see from dependence of S for a ^"^^Pb target that S < 1 when 
< 0.6 MeV. It is thus necessary to determine Sn experimentahy in future in order 
to evaluate Canc and S properly. However, we can say at least that Canc has a 
smaller theoretical error and weaker target dependence than S for any value of B^. 
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Fig. 1. (Color online) dependence of S and Canc- The error bar contains theoretical and 
experimental errors. 



§5. Summary 

We present an accurate method of treating the one-neutron removal reaction at 
intermediate energies induced by both the short-range nuclear and the long-range 
Coulomb interaction. In the theory called ERT, the nuclear and Coulomb breakup 
processes are accurately and consistently treated by CDCC without making the 
adiabatic approximation to the latter, so that the removal cross section never diverges 
and hence reliable even in the presence of the Coulomb interaction. This method is 
applicable also for other inclusive reactions such as the two-nucleon removal reaction 
induced by nuclear and Coulomb interactions. 

ERT is a theory of presenting the Schrodinger equation ()16p to calculate inclusive 
cross sections. When these equations are solved, any accurate method is useful. For 
intermediate and high incident energies, one can use DEA instead of eikonal-CDCC. 
For lower incident energies where the eikonal approximation is not perfectly accurate, 
one should make non-eikonal corrections to inclusive cross sections. This can be done 
easily by using CDCC instead of eikonal-CDCC. 

Canc and S of the last neutron in '^^Ne are evaluated from the one-neutron re- 
moval reaction of the ^^Ne+^^C scattering at 230 MeV/nucleon and the ^^Ne+^'^^Pb 
scattering at 234 MeV/nucleon. Canc has a smaller theoretical error and weaker 
target-dependence than S. Thus, Canc is determined more accurately than S. This 
may change the future strategy for the spectroscopy of unstable nuclei. 



Letters 



9 



When the last neutron in ^^Ne is in the lp3/2 orbit, 5 < 1 for ^ 0.6 MeV, 
and S and Canc have weaker target dependence. When the last neutron is in the 
1/7/2 orbit, meanwhile, 5 > 1 and S and Canc have stronger target dependence. 
These results indicate that the last neutron is mainly in the lp3/2 orbit and imply 
that ^^Ne is deformed. This sort of deformation is expected as an origin of the change 
of magicity.S-EJ 
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